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Abstract. The premetric approach to electrodynamics provides a unified 
description of a wide class of electromagnetic phenomena. In particular, it 
involves axion, dilaton and skewon modifications of the classical electrodynamics. 
This formalism emerges also when the non-minimal coupling between the 
electromagnetic tensor and the torsion of Einstein-Cartan gravity is considered. 
Moreover, the premetric formalism can serve as a general covariant background 
of the electromagnetic properties of anisotropic media. In the current paper, we 
study wave propagation in the premetric electrodynamics. We derive a system of 
characteristic equations corresponded to premetric generalization of the Maxwell 
equation. This singular system is characterized by the adjoint matrix which turns 
to be of a very special form — proportional to a scalar quartic factor. We prove 
that a necessary condition for existance a non-trivial solution of the characteristic 
system is expressed by a unique scalar dispersion relation. In the tangential 
(momentum) space, it determines a fourth order light hypersurface which replaces 
the ordinary light cone of the standard Maxwell theory. We derive an explicit form 
of the covariant dispersion relation and establish it's algebraic and physical origin. 



1. Introduction 

On classical and quantum levels, Maxwell's electrodynamics is a well established 
theory, which results are in a very precise coordination with the experiment. This 
theory, however, can require some principle modifications in order to include non- 
trivial interactions with other physical fields. The following non-complete list indicates 
some directions of the possible alternations: 

• Dilaton field. This scalar partner of the classical electromagnetic field is proposed 
recently as a source of a possible variation of the fine-structure constant [1] , [2] • 

• Axion field. This pseudo-scalar field is believed to play a central role in violation 
of Lorentz and parity invariance [4] — [8j . 

• Birefringence and optical activity of vacuum. These effects are forbidden in the 
standard (minimal coupling) model of interaction between the electromagnetic 
and the gravitational fields. In the Cartan-Einstein model of gravity, the non- 
minimal coupling yields, in general, the nontrivial effects of electromagnetic wave 
propagation [5] — [T2]. 

Although, the mentioned problems belong to rather different branches of the classical 
field theory, their joint treatment can be provided in a unique framework of the 
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premetric electrodynamics. The roots of such an approach can be found in the older 
literature [T^, but its final form was derived only recently, see [11]- [10] and specially 
the book and the references given therein. 

In the premetric construction, the electromagnetic field is considered on a 
bare differential manifold without metric and/or connection. Instead, the manifold 
is assumed to be endowed with a fourth order constitutive pseudo-tensor, which 
provides the constitutive relation between the electromagnetic field strength and 
the electromagnetic excitation tensors. The metric tensor itself is only a secondary 
quantity in this construction. Its explicit form and even its signature is derived from 
the properties of the constitutive tensor, [15] — [21]. 

In the current paper, we study the electromagnetic waves propagation in the 
premetric electrodynamics. From the technical point of view, our approach is similar 
to those used in the relativistic plasmodynamics [21] — [26]. A principle different 
that we are dealing with a metric-free background, thus the norm of the wave covector 
and its scalar product with another covector are not acceptable. Roughly speaking, 
the indices in the tensorial expressions cannot be raised or downed. Although this 
restriction complicates some derivations, it is clearly necessary in the framework 
of the premetric approach. Moreover, it shows that, for the electromagnetic waves 
propagation, the metric tensor is indeed a secondary structure. The metric structure 
must be considered as a result of the properties of the wave propagation and not as a 
predeclared fact. 

The main result of our consideration is a rigorous derivation of the covariant 
dispersion relation. It is shown to be originated in the adjoint matrix of a characteristic 
system of the generalized Maxwell field equations. 

The organization of the paper is as follows: In Section 1, we give a brief account 
of the premetric electrodynamics formalism. Section 2 is devoted to the geometric 
optics approximation and the wave- type ansatz. When this ansatz is substituted 
in the Maxwell system, the former is transformed into a algebraic system of linear 
equations. The algebraic features of this characteristic system is studied in Section 
4. A covariant dispersion relation emerges in Section 5 as a necessary condition for 
existence a non-trivial wave-type solutions in premetric Maxwell electrodynamics. In 
Section 6, we give an application of the developed formalism to the simplest Maxwell 
case. The standard expressions of the Maxwell theory are reinstated. Section 7 is 
devoted to a discussion of the proposed formalism and its possible generalizations. 

2. Premetric electrodynamics formalism 

2.1. Motivations 

In order to represent the motivations of the premetric electrodynamics, we briefly 
recall some electromagnetism models which naturally lead to this generalization. 

A. Vacuum electrodynamics. In the flat Minkowski spacetime, electromagnetic field is 
described by an antisymmetric tensor of the electromagnetic field strength Fij . In an 
orthogonal Cartesian coordinate system {x*} with i — 1,2,3, the dynamics of the 
field is defined from a pair of the first order partial differential equations 

e'^'^^Fjkj =0, F'\j = J\ (2.1) 
Here, the comma denotes the partial differentiation. The Levi-Civita permutation 
pseudotensor e'^'^' with the values {—1,0, 1} is normahzed by e^i^s _ -|^^ 
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The first equation of (|2.ip is completely independent of the metric. In the second 
one, the Minkowski metric, ry*-' = diag(— 1, 1, 1, 1), is involved implicitly. It is used 
here for definition of the covariant components of the field strength, i.e, for raising the 
indices 

F'^ - 77™'7^"i^mn ■ (2.2) 

To have a representation similar to those used below, we rewrite this equation as 

F'^ = ix*^™"F„„, , where x*^'"" = " r?™?7^'™ • (2.3) 

In (j2.ip . the vector field J' describes the electric current. Since the tensor F'-' is 
antisymmetric, the electric charge conservation law 

J\, = (2.4) 

is a straightforward consequence of (|2.ip . 

The relations above are invariant under a subgroup of linear rigid transformations 
of coordinates which preserves the specific form of the Minkowski metric. This group 
includes the instantaneous spatial rotations, Lorentz's transformations and reflections. 



B. Electrodynamics in gravity field. In a non-inertial frame, i.e., in curvilinear 
coordinates on the flat spacetime, the Minkowski metric rf^ is replaced by a generic 
pseudo-Riemannian metric g*-' which components depend on a spacetime point. On 
this background, the transformational requirements are changed. The fleld equations 
must be now invariant under arbitrary smooth transformations of the coordinates. To 
satisfy this transformational requirement, the field equations (|2.ip are modified to 

e'''''F,k,i = , (F^^V^) , = rV~9. (2.5) 



where g = det [gij ) . The covariant components of the electromagnetic fleld strength 
are defined now via a multiplication by the metric tensor components which, in 
contrast to p.2p . depend on a point 

= g'"'g^^Fran . (2.6) 

Observe that now the components of two electromagnetic fields F"^^ and Fmn 
are different functions of a point. In fact they can be treated as two independent 
physical fields. In such an approach, metric tensor comes from a relation between 
these independent fields, i.e., from a physical phenomen. 

Since F^^ is antisymmetric, the inhomogeneous field equation of (12. 5|) yields a 
modified electric charge conservation law 

(J^V^),,, = 0. (2.7) 

In fact, this equation is not a conservation law for the vector field ,P itself. What is 
really conserved is the product of J* with the root of the metric determinant. It means 
that a conserved electric current cannot be described by a covariant vector field so a 
redefinition of this basic notion is necessary. Instead of treating it as a vector field, 
the electric current has to be considered as a weight (+1) pseudo-vector density field 

T = J'^f~g. (2.8) 

Also a weight (+1) pseudo-tensor density of electromagnetic excitation 

= F'^y^ (2.9) 
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has to be involved. Under smooth transformations of the coordinates — > with 
the Jacobian L = detidx'^ /dx^), the transformation law for these pseudo-tensorial 
quantities involves an additional factor In order to have a covariant field 

equation, this factor must be compensated. The first order partial derivatives of 
the term \/—g makes the job and the whole equation is covariant. 

Consequently, the general covariant field equations take the form 

e'^'^'i^jM =0, W\j=J\ (2.10) 
while the general covariant charge conservation law is written as 

J\^^0. (2.11) 
The constitutive relation between two basic fields takes the form 

F'^ = ^x''"'"F„,n (2.12) 

where the constitutive pseudotensor 

X =V-9[9 9 -9 9 ) (2-13) 

is involved. 

Although the described modification serves the curvilinear coordinates on a flat 
manifold, it is well known to be enough also for description of the electromagnetic 
field in a curved spacetime of GR. Both field equations (|2.10p and the conservation 
law (|2.1ip are general covariant even being written via the ordinary partial derivatives. 



C. Electrodynamics in anisotropic media. For an anisotropic media in a flat 
Minkowski space, Maxwell's electrodynamics is described by two pairs of 3D vectors 
Ea,B°' and Z?",i?Q, where the Greek indices are assumed to obtain the spatial 
values, a,/3, • • • = 1,2,3. In the 4D notation, these vectors are assembled into two 
antisymmetric tensors: the electromagnetic strength tensor Fij with the components 

Foa — Ea , Fap — ~€af3jB'^ , (2.14) 

and the electromagnetic excitation tensor H^^ with the components 

ij"" = , H"'^ = e"f^''H^ . (2.15) 

In the 4D notation, the Maxwell field equations for the electromagnetic field in an 
anisotropic media are written in the form 

e'^'^'i^^j.fc =0, H'\,^J\ (2.16) 

An additional ingredient, the constitutive relation between two electromagnetic 
tensors, Fij and if'-', describes the characteristic properties of the media. For a 
wide range of anisotropic materials, a linear constitutive relation is a sufficiently good 
approximation 

^e'^^'Ep + ^'^pB" , Ha=li-lBf' +^jEp. (2.17) 

The electromagnetic current conservation law = is a consequence of the field 
equation (|2.16p . Notice also that the equations (|2.16p are invariant under arbitrary 
constant linear transformations of the coordinates. 
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2.2. Premetric field equations 

The models accounted above show some similarity: 

• The electromagnetic field is described by two 2-nd order antisymmetric tensors. 
In differential form notation, the field is represented by two 2-nd order differential 
forms — one twisted and one untwisted. 

• All the models are described by similar systems of two first order partial 
differential equations. 

• The antisymmetric tensorial fields are related by a linear constitutive relation. 

• Even in vacuum electrodynamics, the metric of the manifold emerges only via a 
special 4-component tensor, i.e., it plays only a secondary role. 

The accounted similarity naturally leads to a premetric generalization of the classical 
electrodynamics. For a comprehensive account of this subject, see [19], [20] and the 
references given therein. 

In the premetric approach, two differential field equations for two second 
order differential forms, the electromagnetic field strength F and the electromagnetic 
excitation Ti are postulated: 

dF = 0, dn = J. (2.18) 

In (|2.18|) . F is an even (untwisted) differential form. It does not change under arbitrary 
transformations of coordinates. Alternatively, Ti and J are odd (twisted) differential 
forms. Under a change of coordinates with a Jacobian L = det {L^j), they multiplied 
by the sign factor of L. Namely such identification of the electromagnetic fields 
guarantees the proper integral conservation laws for magnetic flux and electric current, 
see [19] . Notice that the domain of integration for a twisted form must carry a chosen 
orientation. 

Both equations (|2.18p are expressed via differential forms thus they are manifestly 
invariant under arbitrary smooth coordinate transformations. In a more general 
setting, see [19], these equations can be considered as consequences of two integral 
conservation laws: one for the magnetic ffux and one for the electric current. 

In a coordinate chart, we represent the forms as 

F=^ F,jdx' /\dx' , H = i H*^ey„„ da;" A dx" , (2.19) 

while 

J = 1 J'ey„„ dx^ A dx™ A dx" . (2.20) 

Thus, the components Fij constitute an ordinary antisymmetric tensor while the 
components and JT"* are pseudotensor densities of weight (-1-1). 

Applying the exterior derivatives to (|2.19p we rewrite the field equations (|2.18p 
in the tensorial form 

e'^'^'Fjkj^O, W\j=J\ (2.21) 

Even being written via the ordinary partial derivatives, these equations are covariant 
under arbitrary smooth transformations of coordinates. Notice also that a metric 
tensor or a connection are not involved in the construction above. One can even say 
that these structures are not defined (yet) on the space. In this sense, the construction 
is premetric. Particularly, in such approach, the covariant components of the field 
strength tensor and the contravariant components of the excitation tensor cannot be 
introduced — the indices cannot being raised or downed. 
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2.3. Constitutive relation 

The system (|2.2ip involves 8 equations for 12 independent variables so it is 
undetermined. Moreover, the 2-form H describes a field generated by a charged 
source, while the 2-form F describes some other field which acts on a test charge. 
On a current stage of the construction, these two fields are formal and completely 
independent. It means that an interaction between two charges is not yet involved. 
In order to close the system and to involve an interaction, a constitutive relation 
between the fields Fij and TV'^ must be implicated. The simplest choice of a local 
linear homogeneous relation 

= \x'''''Fu (2.22) 

is wide enough to describe the most observation data of the ordinary electrodynamics 
and even involves some additional electromagnetic effects, such as axion, dilaton and 
skewon partners of photon, see [12]. Also for the electromagnetism into the non- 
magnetized media, the linear constitutive relation is a good approximation. For the 
nonlinear extensions of the premetric approach, see |23| . The non-local constitutive 
relations was considered recently in [55] . 

Recall that the physical space is considered as a bare manifold without metric 
or connection. All the information on its geometry is encoded into the constitutive 
pseudotensor x^''^^ which can depend of the time and position coordinates. By 
definition, this pseudotensor inherits the symmetries of the antisymmetric tensors 
Fij, T-C^ . In particular, 

yAjkl _ ^[ij]kl _ ^ij[kl] ^ (2.23) 

Thus, in general, the fourth order constitutive tensor x*-''^' has 36 independent 
components instead of 4"^. Due to the Young diagrams analysis, under the group of 
linear transformations, such a tensor irreducible decomposed into three independent 
pieces 

^ijkl _ (l)^ijkl {2)^ijkl {3)^ijkl ^ (2.24) 

The axion part (1 component) and the skewon part (15 components) are defined 
respectively as 

(3)^ijkl ^ ^tjkl] ^ (2)^ijkl ^ 1 (x*-'''' - X*"''^) • (2.25) 

The remainder (^)x*"''^' is a principal part of 20 independent components. One can also 
extract a scalar factor from the principle part which was identified recently ^22, with 
the dilaton partner of electromagnetic field [T] , [5] . 



2.4^. Uniqueness and existence problem 

In this paper, we are studying the system of field equations (|2.2ip with a linear 
constitutive relation l|2.22p . The constitutive tensor x^^^^ is irreducible decomposed 
to a sum of three independent pieces. In absence of a metric tensor, a successive 
decomposition is not available. Since the decomposition (|2.24p is irreducible, the 
difference pieces of the constitutive tensor can appear alone or in certain linear 
combinations. Thus, in principle, we have several non-trivial electromagnetic models. 
Not all of them, however, confirm with the physical meaning of the field equations 
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In order to extract some physics meaningful models, we will apply the following 
requirements: 

(i) Hyperbolicity. The free field must propagate in spacetime by waves. Thus 
the partial differential equations has to be of the hyperbolic type. 

(ii) Existence and uniqueness. For physically meaningful sources, the system 
(|2.2ip with the relation (|2.22p can serve as a physical model only if it has a solution 
and this solution is unique. 

It is easy to see that these requirements are not trivial. In particular, (ii) is 
not satisfied in the following case: Let only the axion part ^^'>^^i'^^ be nonzero. Let 
it also be a constant. In this case, the equations (|2.21[) with the relation (|2.22[) are 
inconsistent for J^' 7^ and undetermined for J'^ — 0. Consequently an isotropic 
and homogeneous pure axionic media is forbidden by the uniqueness and existence 
conditions. 

One of the main tasks of the current paper is to establish under what conditions 
a generic constitutive tensor can represent a meaningful physical system. 

3. Approximation and ansatz 
3.1. Semi- covariant approximation 

When the constitutive relation (|2.22p is substituted into the field equation (|2.21l, ). we 
remain with 

lx'"''Fki., + lx'''\,Fu=r. (3.1) 

The first term here describes how the electromagnetic field changes in a spacetime of 
constant media characteristics. Alternatively, the second term describes the temp of 
the spacetime variation of the media characteristics for a constant electromagnetic 
field. In this paper, we restrict to the geometrical approximation of optics. In 
particular, we neglect with the second term of (|3.1[) relative to the first one. In other 
words, we restrict to a media which characteristics change slowly on the characteristic 
distances of the change of the fields. Note that such approximation is not always 
applicable. In particular, the Carroll-Field- Jackiw modification of electrodynamics 
^\ , can be reformulated as a premetric electrodynamics [3] , [7] , [5] where the first term 
of (|3.ip vanishes, while the birefringence effect comes namely from the second term. 

In the framework of the geometrical approximation, we remain with a system of 
8 equations for 6 independent components of the electromagnetic field strength Fij 

In a special case of the Maxwell constitutive tensor, the approximation used here yields 
the inhomogeneous field equation of the form 

g'^g'^Fki,, = r ■ (3.3) 

This is an approximation of the covariant equation (|2.5p when the derivatives of 
the metric tensor are considered to be small relative to the derivatives of the 
electromagnetic field. Such a semi- covariant approximation is preserved for arbitrary 
coordinate transformations with small spacetime derivatives. 
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3.2. Eikonal ansatz 

To describe the wave- type solutions of the field equations (j3.2p . we consider an eikonal 
ansatz. Let the electric current be given in the form 

r^fe\ (3.4) 

and let the corresponded field strength be expressed as 

F^j = /ye" . (3.5) 

Here the eikonal cr is a scalar function of a spacetimc point. The tensors and 
fij are assumed to be slow functions of a point. The derivatives of a give the main 
contributions to the field equations. Define the wave covector 

In this approximation, the conservation law for the electric current J'_,; ~ takes a 
form of an algebraic equation 

rm = 0. (3.7) 
Substituting (|3.41 13. 5p into the field equations (|3.2p and removing the derivatives of the 
amplitudes relative to the derivative of the eikonal function, we come to an algebraic 
system 

e^'^'q^fki^O, x'"''q,fM=1f ■ (3.8) 
The same system was derived in [19] by mean of Hadamar's discontinuity propagation 
method. This fact indicates that a simple approximation used here is not less general 
than the one used in [19] . 

Observe a remarkable property of (|3.8p . If all the quantities involved here 
are assumed to transform by ordinary tensorial transformation rules with arbitrary 
pointwise matrices, both equations are preserved. In other word, these equations are 
straightforward expanded to a general covariant system. This is in spite of the fact 
that the approximations used in their derivation are not covariant. This property 
is generic for quasilinear systems which leading terms (the higher order derivatives 
expressions) are linear and thus preserve their form even under arbitrary pointwise 
transformations . 



4. Characteristic equations 

^.1. The linear system 

The approximation (|3.2p and the wave-type ansatz (|3.41 13. 5p yield a linear system 
p.8p of 8 equations for 6 independent variables. This algebraic system will serve as 
a starting point of our analysis. Observe first that the system is not overdetermined. 
Indeed, when both equations are multiplied by a covector g^, they turn to trivial 
identities, provided that the electric current is conserved. Thus we have two linear 
relations between 8 linear equations p.Sp for 6 independent variables. It means that 
the rank of the system p.8p is less or equal to six. The physical meaning of this system 
requires the rank to be exactly equal to 6. Indeed, the unknowns fki of this system are 
physically measurable quantities. Thus, for an arbitrary conserved current j*, they 
have to be determined from ()3.8p uniquely. This physical requirement puts a strong 
algebraic constrain on the system (|3.8p — it coefficients must form a matrix of a rank 
of 6. In fact, it is a constrain on the components of the constitutive pseudotensor 
^ijki ^ formal expression we will derive in sequel. 
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4-2. The homogeneous equation 

The homogeneous equation (|3.8p is exactly the same as in the standard MaxweU 
theory. We give here a precise treatment of this equation mostly in order to establish 
the notation and to illustrate the method used in sequel. 

Proposition 1. A most general solution of a linear system 

e'^^'q^fki = (4.1) 

is expressed as 

fkl = ^ {akQi - aiQk) , (4.2) 
where au is an arbitrary covector. 

Proof: The expression (|4.2p is evidently a solution of (|4.ip . In order to prove that 
it is a most general solution, we first notice that (|4.ip is a linear system of 4 equations 
for 6 independent variables fij. The 4x6 matrix of this system a*'^' = e^^^^qj with 
ik = 01, 02, 03, 12, 23, 31 and ? = 0, 1, 2, 3 is given by 



( ' 
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(4.3) 



The rows of this matrix satisfy a linear relation 

a'^^'g, = 0. (4.4) 

Thus its rank is of 3 or less. If an arbitrary row is removed now from (|4.3p . the 
remaining three columns are assembled in the echelon form. Thus the matrix (j4.3p 
has exactly a rank of 3. Consequently, a general solution of (14. ip has to involve 
6 — 3 = 3 independent parameters. This is exactly what is given in (14. 2p . Indeed, 
although the arbitrary covector Oi has 4 independent components, only 3 of them 
are involved in (|4.2p . In particular, the vector Oi proportional to qi does not give a 
contribution. Thus, ()4.2p is a most general covariant solution of (|4.ip . ■ 

It is well known that the homogeneous field equation (|3.2p is solved in term of 
the standard vector potential Ai. The covector appeared in (|4.2p is similar to the 
Fourier transform of of the potential Ai. As ususal, this covector is arbirtrary if only 
the homogeneous field equation is taken into account. 



4.3. The inhomogeneous equation 

Let us turn now to the inhomogeneous equation of (|3.8p . Substituting the solution 
(|4.2p . we arrive to an algebraic system 

x'^'^rnqjak^f. (4.5) 

Observe that this is a system of 4 equations for 4 variables a^. The matrix of this 
system 

M''= = x'^'^'qiq, , (4.6) 

will be refereed to as a characteristic matrix. We will see in sequel that the wave 
propagation depends exactly on the specific combination of the components of the 
constitutive pseudotensor x*^*^' which are involved in M**^. 



On light propagation in premetric electrodynamics. Covariant dispersion relation 10 



When the irreducible decomposition (|2.24p is substituted into the characteristic 
matrix M'^ , the completely antisymmetric axion part (3);^'^'=' evidently does not 
contribute. As for the other two pieces, the principle part is involved only in the 
symmetric part of the matrix M*-', while the skewon part is involved only in its 
antisymmetric part. Formally, we can write 

So, the characteristic matrix M^^ is irreducibly decomposed as 

M'^ = Af(*'=^) ((i)x) + A^''''' (^^^x) • (4.8) 
In the characteristic matrix notation, equation (j4.5p takes the form 

M^'^Ofc = . (4.9) 

The following two facts will play an important role in our analysis: 

(1) Gauge invariant condition. Due to the antisymmetry of the constitutive 
pseudotensor x^''^^ in its last two indices, an identity 

M'^qk = (4.10) 

holds true. It is a linear relation between the rows of the matrix AP^ . It means that 
every solution of (|4.5|) is defined only up to an addition of a term a.; ~ qi. This addition 
is evidently unphysical since it does not contribute to the electromagnetic strength. 
Consequently, relation (j4.10p has to be interpreted as a gauge invariant condition. 

(2) Charge conservation condition. Another evident identity for the matrix M**^ 
emerges from the antisymmetry of the constitutive pseudotensor ^^^^^ in its first two 
indices 

M'^q, = Q. (4.11) 

It is a linear relations between the columns of the matrix M^^ . Being compared with 
(|4.9p . it yields j'^qi = 0. Consequently, the relation (|4.1ip has to be interpreted as a 
charge conservation condition. 

Thus we arrive to a some type of a duality between the charge conservation and 
the gauge invariance. Note that this duality is expressed by a standard algebraic fact: 
For any matrix, the column rank and the row rank are equal one to another. 

Due to the conditions indicated above, the rows (and the columns) of the matrix 
M'^ are linearly dependent, so its determinant is equal to zero. It can be checked 
straightforwardly, but one has to apply here rather tedious calculations. 

5. Dispersion relation 

5.1. How it emerges 

In the vacuum case of the free electromagnetic waves, (j4.5l) takes a form of a linear 
homogeneous system of four equations for four components of the covector a^, 

X^'^^'qiq^au^Q, or M^'^a,- = . (5.1) 

The gauge relation (|4.10[) can be interpreted as a fact that 

ai = Cqi (5.2) 

is a formal solution of (|5.ip . This solution does not contribute to the electromagnetic 
field strength so it is unphysical. Hence, the formal system (|5.ip can describe a 
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nontrivial wave propagation, only if it has an additional solution which must be linearly 
independent on (|5.2p . Consequently (|5.ip must have at least two linear independent 
solutions. A known fact from the linear algebra that a linear system has two (or more) 
linearly independent solutions if and only if the rank of the matrix M*-' is of two (or 
less). In this case, the adjoint matrix (constructed from the cofactors of M*-') is equal 
to zero, Aij = 0. 

In order to present a formal expression of this fact we will start with a formula 
for the determinant of an arbitrary 4-th order matrix, 

det{M) = ie,,,,,,3e,,,,,,3M*^M^^^'W^^-^M'«"^ . (5.3) 

Here the standard convention £0123 = 1 is accepted. The components of the adjoint 
matrix are expressed by the derivatives of the determinant relative to the entries of 
the matrix, 

_ddet{Ml 

dM^i ■ ^^-^^ 

Explicitly, 

1 

3!' 

Consequently, we derived a physically motivated condition on the components of the 
constitutive pseudotensor Xijki- 



A^, = T^^^^^^.^,e,,,,,nM''''M^^^-M^-^- . (5.5) 



Theorem 2. The Maxwell system with a general linear constitutive relation has 
a non-trivial wave-type solution if and only if the adjoint of the matrix M^^ = x^-'^^qiqj 
is equal to zero, i.e., 

A,, - , (5.6) 

or explicitly, 

e»ii.i3e«i...3M^"'^A^*^^'^M"-''^ - . (5.7) 

Since the adjoint matrix has 16 independent components, it seems that we have to 
require, in general, 16 independent conditions. In fact, the situation is rather simpler. 
The following algebraic fact is important in our analysis, so we present here its formal 
proof. 

Proposition 3. Let a square n x n matrix NP^ satisfies the relations 

lVP^q, = Q, A'r^qj=0, (5.8) 

for an arbitrary nonzero n-covector qi. The adjoint matrix Aij — Adj{KP^) is 
proportional to the tensor square of qi, i.e., 

Aij = \{q)qiqj . (5.9) 

Proof: Due to (|5.8p . the rows (and the columns) of AP^ are linear dependent, 
so the matrix is singular and its rank is equal to n — 1 or less. If the rank is less than 
(n — 1), the adjoint matrix is identically zero and (|5.9p is satisfied trivially for A = 0. 

Let the rank of M^^ be equal to (n — 1). In this case, qi is a unique covector 
(up to a multiplication on a constant), that satisfies (|5.8p . It is well known that for a 
matrix of a rank of (n — 1), the adjoint Aij is a matrix of a rank of one. Moreover, an 
arbitrary rank one matrix can be written as a tensor product of two covectors 

Aij = UiVj . (5.10) 
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Let us show now that both these covectors must be proportional to qi. Indeed, the 
product of an arbitrary matrix with its adjoint is equal to the determinant of the 
matrix times the unit matrix (this is the generalized Laplace expansion theorem). In 
our case, M*-' is a singular matrix so 

Aij M''^ = A,j m''^ = 0. (5.11) 

Substituting here (|5.10p . we have the relations 

UiVjM'^ = u^VjM^'' = , (5.12) 

or, equivalently, 

UrM'^ = VjM^'' = 0. (5.13) 

Comparing this pair of relations to (j5.9p and remembering that qi is unique up 
to a multiplication on a constant, we conclude that both covectors Ui and Vi are 
proportional to qi. Thus (|5.4[) indeed takes the form (|5.9p . ■ 

Consequently, in order to have a physically non-trivial vacuum wave-type solution, 
the system (15. ip has to satisfy a unique scalar condition 

X{q)^0. (5.14) 

In fact, this is an expression for the principle dispersion relation. 



5. 2. Some basic properties of the dispersion relation 

Even without an explicit expression for the function A(q), we are ready now to 
derive some characteristic properties of the dispersion relation (|5.14p . Some of these 
properties was recently derived in [19] by involved straighforward calculations. In 
our approach, these properties are immediate consequences of the definition of the 
A-function (jO)) . 

Corollary 4. X{q) is a homogeneous 4-th order polynomial of the wave covector 
Qi, i-e., 

m = g'"''q.mkqi , (5.15) 

where Q^^^^ is a pseudotensor independent on qi . 

Indeed, the adjoint matrix is a homogeneous polynomial of the six order in qi. By 
(|5.9p . after extracting the product g^gj we remain with a sum of terms every one of 
which is a product of four components of the covector g^. Since A(g) is a pseudoscalar, 
gijki ^ pseudotensor. 

Corollary 5. A(q) is a homogeneous 3-rd order polynomial of the constitutive 
pseudotensor ^^i'^^ i.e. it is of the form 

g^'"' = [ee ■ ixXX)]''"' , (5.16) 
where the dot symbol denotes the tensor contraction with the Levi-Civita pseudotensor 
e. 

Indeed, the adjoint matrix is a sum of terms which are cubic in the matrix M*^ . 
Every such a term is a product of three x's which remains also after extracting the 
product qiqj in the right hand side of l|5.9p . 

Corollary 6. The equation X{q) — defines a two-folded algebraic cone. 
For a given set of vectors K — {xi,X2,- ■ ■}, an algebraic cone is defined as a 
set of vectors CK = {Cxi,Cx2, ■ ■ ■}, where C > is an arbitrary number. Due to 
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the homogeneity of the dispersion relation (|5.14p . every its solution is defined up to 
a product on a constant. Thus we have two folds of the cone — one for C > and 
the second for C < 0. This two-folded algebraic cone is a prototype of a light cone 
emergered by an Lorentz metric of vacuum electrodynamics or by optical metric of 
electromagnetism in a dielectric media. 

Corollary 7. The axion part of the constitutive tensor does not contribute to the 
function X{q). In other words, 

a(«X+(^^X+<^'x)-a(Wx+(^)x) . (5.17) 

Indeed, the axion part does not contribute to the matrix AP^ , so it does not 
appear in its adjoint. 

Corollary 8. The skewon part alone does not emerges a non-trivial dispersion 
relation. In other words, a relation 

X(^'\)^0 (5.18) 

holds identically. 

Indeed, in order to have a non-trivial (non-zero) expression for A(g), the rank 
of the matrix M'^ has to be equal to three. The skewon par alone generates an 
antisymmetric matrix Since the rank of an arbitrary antisymmetric matrix is 

even, the skewon part alone does not emerges a non-trivial dispersion relation. 

Corollary 9. A non-trivial (non-zero) dispersion relation emerges only if the 
principle part of the constitutive tensor is non-zero 

(i)xy^O. (5.19) 

This is an immidiate result of the previous statements. 



6. Dispersion relation in an explicit form 
6.1. Covariant dispersion relation I 

Our task now is to derive an explicit expression for the dispersion relation. Recall 
that it is represented by a scalar equation 

A(g)=0, (6.1) 

where the function X{q) satisfies the equation 

Aij = \{q)qiqj (6.2) 

for the adjoint matrix Aij of the characteristic matrix M^^ . To have an explicit 
expression for X{q), it is necessary "to divide" both sides of (|6.2p by the product qiqj. 
Certainly such a "division" must be produced in a covariant way. We will look first 
for a solution of this problem in a special coordinate system. Let a zeroth (time) axis 
be directed as a wave covector, i.e., go = 9j <?! = 92 = 93 = 0. Substituting into (|6.2p 
we have 

Hq)q' - i;eo^,^,^seonJ.JsM^'''M^^^^^p■■^^^ . (6.3) 
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Due to the symmetry properties of the Levi-Civita pseudotensor, aU four-dimensional 
indices can be replaced by the three-dimensional ones (a, /3 = 1,2,3). So we get 

KiW = ^eoaia.aaeoft&fe^"''''^"''''^"'* ■ (6.4) 
In the chosen system, the non-zero components of the matrix NP^ are 

M<-P = x""''"g™gn = X""^"(7' . (6.5) 

Consequently, 

1 
3! 

where the three dimensional Levi-Civita pseudotensor eojojag = coaioaaa is involved. 
The non-covariant dispersion relation takes the form 19] 

^o.a.aa^^.fcfeX"^"''^ V"^°''^°X"^°''^°9' = . (6.7) 

Note that this expression is only formal. Indeed, it is equivalent to go = 9 = Oj i-e, 
all the components of the wave covector are equal to zero. This statement cannot 
represent any physical meaningful property of the wave propagation. It is due to 
the fact that the system with a zero wave covector cannot be realized by any matter 
particles, i.e., there is not an observer with a 4- velocity equal to the wave covector. 

In [TS], [12], the equation (|6.7p was derived by consideration of the three- 
dimensional determinant of the system. In spite of the problems with it's physical 
neaning, the expression (|6.7[) was used as a guiding form for a proper dispersion 
relation. In particular, by physical intuition, this equation was generalized to a 
covariant four-dimensional dispersion relation 
1 
4! 

The A-function can be read off from this equation as 

m - ^e.n...3£j.u2.3x"^^'''x''*^^'^^x''^^'^^'^9agbgc<zd • (6.9) 

Although, from a formal point of view, the described way of the derivation of the 
covariant dispersion relation is very much doubtful, the final result (|6.8p turns to be 
correct, as we will show in sequel. 

6.2. Covariant dispersion relation II 

We will give now a pure covariant derivation of the scalar function \{q) and of the 
corresponded dispersion relation. Differentiation of (|5.9[) relative to the components 
of the covector yields 



^u^^,^,^,n3.3.X'"'"'x'"''''x''''''''qaqbqcqd = . (6.8) 



dA,j _ dX{q) 



q,q,+\[q){5Tqj+5fq,). (6.10) 



Let us contract this equation over the indices m and i and use the Euler formula for 
a fourth order homogeneous function X{q). Consequently, we derive 



<d\{q)qj . (6.11) 
ression is given by 

q, + msA . (6.12) 



dqi 

A second order derivative of this expression is given by 
cfA,, ^^(d\{q) 



dqidqm \ dq 
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Summing now over the indices m and j and using once more the Euler formula, we 
derive 

Consequently, we have proved the following 

Theorem 10. For the Maxwell system with a general local linear constitutive 
relation, the dispersion relation is given by 

rP- A 

1^ = 0. (6.14) 

In order to have an expression of the A-function in terms of the matrix M*-' , we 
calculate the derivatives of the adjoint matrix 

Substituting into (|6.14p and calculating the second order derivative, we get 

A(.) = ^e.....3e,.,.,3 (^^M--^ + '^^) • ^'-''^ 

This expression may be useful for actual calculations of the dispersion relation for 
different media. 

An explicit expression of the A-function via the constitutive pseudotensor is 
calculated by the derivatives of the matrix 

M'^ = X^'^Qaqb = -X"""qaqb - -x'^"'^' QaQb ■ (6.17) 
The first order derivative is given by 

«)\/filil f) 

-7^ = 7^ (X'^"^'^"9m9n) - -2x''^^^^''qm . (6.18) 

Hence, the second order derivative reads 

-jTir- = -2x''^''^'' ■ (6.19) 
Consequently, the left hand side of (|6.16p takes the form 
A(g) = ^eu,.,^,e,,,,,,,{-x''^''^''M^^^^+Ax'^^^^^^^ 



3-4! 

— ^ r ■ f - ■ ■ f ^i2{ab)j2 I A ii{ia)ji i2{jb)j2\ i3{cd)j3 „, „ „ , 

— 3 41 *i*2*3 '^JJlj2j3 A < A J A HaHbHc'id ■ 

(6.20) 

We finally have the covariant dispersion relation in the form 

^^^^^2^^^n^3.n (x'^^^^'^'' X'^"'"' + ^X''^'^^'' X'''^"''^'')x'''^'' ^aqbqcqd = . (6.21) 

This equation turns out to be equivalent to (|6.8p . The direct proof of this fact was 
provided by Yu. Obukhov by a sequence of rather involved algebraic manipulations 



On light propagation in premetric electrodynamics. Covariant dispersion relation 16 

The following decomposition represents the contribution of the skewon part into 
the dispersion relation. Different forms of it can be found in jlS] . 

Proposition 12. Due to the irreducible decomposition of the constitutive pseudo- 
tensor, the dispersion relation \{x) = is given by 

1 

' dqidqj 



A (wx) + r---^--- [^"'""^ i'"'^) {'"'^) ('''^) 



= 
(6.22) 



Proof: The relation is follows straightforwardly when the decomposition (|4.8[) is 
substituted into (|6.14p and the antisymmetric terms are removed. I 

7. MELXwell electrodynamics reinstated 

7.1. Maxwell constitutive pseudotensor 

When the premetric scheme is applied on a manifold with a prescribed metric tensor 
gij, the standard Maxwell electrodynamics (|2.5p is reinstated by substitutions 

= x/^.g""5^"^^™„ ^IV^ - ^mn (7. 1) 



2 

Recall that the electric current vector density is expressed as 

J'^r^, (7.2) 

The constitutive relation (|7.ip corresponds to a choice of a special Maxwell-Lorentz 
constitutive pseudotensor 

We use here a system of units in which a constant with the dimension of an admittance 
(denoted by Aq in [l^) is taken to be equal to one. When (|7.3p is substituted in (|2.2ip 
we return to the standard Maxwell electrodynamics system (|2.5p . 

7.2. Dispersion relation 

Let us derive the standard metrical dispersion relation. The characteristic matrix 
corresponded to the constitutive pseudotensor (|7.3p takes the form 

M^' = V^(5V-9^9^) , (7.4) 

where the notations — g^-'qiqj and q^ = g^"^qm are used. The adjoint of this matrix 
is calculated straightforwardly 

= l^ VW^^n^.^s^m,.,. {9'''' l' ~ 1''^'') C?^^'^?' - q'^ q'^) {g'^'^q' ~ ^t'q'') 
= l^ VW^^n^.^s^m,.,. {9'''' 9'^'^ 9'''' q' " 3g'^'^ g'^'^ q^W' q') 

[-'Q9tjq^ + Hgij9iiji - gi3i9iij)q'"q^'q'*^] ■ (7-5) 



3! 

The first two terms are cancel one another so the adjoint matrix remains in the form 

= -VWq^q3q'^ ■ (7-6) 
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(7.7) 



and the dispersion relation takes the standard form 



,2 







g'^^q^qj = 0. 



(7.8) 



From (|7.7p we also deduce that the pseudotensor appeared in (|5.15l5.16p takes the 
form 



8. Results and discussion 

Premetric electrodynamics can be viewed as a general framework for description of a 
wide class of electromagnetic effects. In this paper we discussed the geometric optics 
approximation of the wave propagation in this model. Due to the standard procedure 
of partial differential equations theory, such approximation represents the leading 
contribution to the corresponded solutions. We derived a covariant dipersion relation 
and showed that this relation represents the existance of the wave-type solution in 
the premetric electromagnetic model. It should be noted that our expression of 
the covariant dispersion relation is not less complicated that the one represented in 
literature, [19] . An advantage of our approach is that we give an straightforward 
covariant procedure how the dispersion relation can be derived for various constitutive 
relations. For this, one does not need to deal with the explicit covariant formula at all. 
It is enough to construct the characteristic matrix M*-' and calculate it's adjoint A^j. 
Due to the proposition, proved in the paper, the extra factors qiqj are separated from 
Aij and the remain part is the essential term of the dispersion relation. We have shown, 
how this procedure works in the case of the standard Maxwell constitutive relation. 
The problem of uniqueness of the wave-type solution in the premetric electromagnetic 
model is related to another principle notion — the photon propagator, [28j, [8]. A 
detailed consideration of this quantity and of it's relation to the uniqueness problem 
will be represented in a separate publication. 

Acknow^ledgment 

I would like to thank Friedrich Hehl and Yuri Obukhov for most fruitful discussions. 
References 

[1] J. D. Bekenstein, Phys. Rev. D 25 (1982) 1527. 

[2] J. D. Bekenstein, Phys. Rev. D 66, 123514 (2002) 

[3] F. Wilczek, Phys. Rev. Lett. 58, 1799 (1987). 

[4] S. M. Carroll, G. B. Field and R. Jackiw, Phys. Rev. D 41, 1231 (1990). 

[5] V. A. Kostelecky R. Lehnert and M. J. Perry Phys. Rev. D 68, 123511 (2003) 

[6] Y. Itin, Phys. Rev. D 70, 025012 (2004) 

[7] Y. Itin, Gen. Rel. Grav. 40, 1219 (2008) 

[8] Y. Itin, Phys. Rev. D 76, 087505 (2007) 

[9] O. Preuss, M. P. Haugan, S. K. Solanki and S. Jordan, Phys. Rev. D 70, 067101 (2004) 

[10] G. F. Rubilar, Y. N. Obukhov and F. W. Hehl, Class. Quant. Grav. 20, L185 (2003) 

[11] Y. Itin and F. W. Hehl, Phys. Rev. D 68, 127701 (2003) 

[12] O. Preuss, S. K. Solanki, M. P. Haugan and S. Jordan, Phys. Rev. D 72, 042001 (2005) 




(7.9) 



On light propagation in premetric electrodynamics. Covariant dispersion relation 18 

[13] E.J. Post, Formal Structure of Electromagnetics - General Covariance and Electromagnetics 

(North Holland, Amsterdam, 1962, and Dover, Mineola, New York, 1997). 
[14] G.F. Rubilar, Yu.N. Obukhov, and F.W. Hehl, Int. J. Mod. Phys. Dll, 1227 (2002) 
[15] Y. N. Obukhov, T. Pukui and G. F. Rubilar, Phys. Rev. D 62, 044050 (2000). 
[16] C. Lammerzahl and F. W. Hehl, Phys. Rev. D 70, 105022 (2004). 
[17] F. W. Hehl and Y. N. Obukhov, Lcct. Notes Phys. 702, 163 (2006) 
[18] Y. Itin, Phys. Rev. D 72, 087502 (2005) 

[19] F.W. Hehl and Yu.N. Obukhov, Foundations of Classical Electrodynamics: Charge, Flux, and 

Metric (Birkhiiuscr, Boston, MA, 2003). 
[20] F.W. Hehl, Y. Itin, Yu.N. Obukhov, Recent developments in premetric classical electrodynamics, 

in Notebooks on Physical Sciences XVIII: Conferences, Al (2005) 375-408 (Institute of 

Physics: Belgrade, 2005) [arXiv.org/physics/0610221]. 
[21] Y. Itin and F. W. Hehl, Annals Phys. (N.Y.) 312, 60 (2004). 
[22] F. W. Hehl and Yu. N. Obukhov, Gen. Rel. Grav. 40, 1239 (2008) 
[23] Y. N. Obukhov and G. F. Rubilar, Phys. Rev. D 66, 024042 (2002) 
[24] D. B. Melrose, Plasma Phys. 15, 99 (1973) 
[25] M. Gedalin, D. B. Mehose, Phys. Rev. E 64, 027401 (2001) 

[26] D. B. Melrose, Quantum plasmadynamics: unmagnetized plasmas, Lect. Notes Phys. 735 (2008) 
[27] Y. N. Obukhov. Private communication. 
[28] Y. Itin, J. Phys. A 40, F737 (2007) 



